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ABSTRACT 

We have used high-resolution spectroscopy to observe the Kepler-16 eclipsing binary as a double- 
lined system, and measure precise radial velocities for both stellar components. These velocities yield a 
dynamical mass-ratio of (j = 0.2994±0.0031. Wh en combined with t he inclination, i = 90?340lto °oi9, 
measured from the Kepler photometric data by iDovle et alj ()2011[ ). we derive dynamical masses for 
the Kepler-16 components of Ma = 0.654 ± 0.017 Mq and Mb = 0.1959 ± 0.0031 Mq, a precision of 
2.5% and 1.5% respectively. Our results confirm at the ~2% level the mass-ratio derived by Dll 
with their photometric-dynamical model, q = 0.2937 ± 0.0006. These are among the most precise 
spectroscopic dynamical masses ever measured for low-mass stars, and provide an important direct 
test of the results from the photometric-dynamical modeling technique. 

Subject headings: binaries: eclipsing — binaries: spectroscopic — stars: fundamental parameters — 
stars: individual (Kepler-16) — stars: low-mass — techniques: radial velocities 



1. INTRODUCTION 

The NASA Kepler space mission is monitoring about 
150,000 stars with a photometric precision of a few parts 
per mill ion in an effort to find transiting Earth mass 
planets (|Borucki et al.|[2010D . This program has lead to 
the discovery of more than 2000 eclipsing binaries (EBs) 
in t he Kepler field, for which public light- curves now ex- 
ist (|Prsa et al.l 120111 : iSlawson et al.ll2011h . These light- 
curves constrain the inclination angle, i, of the orbit, 
and when coupled with radial velocities (RVs) derived 
from double-lined spectroscopic observations (SB2), they 
yield precise dynamical masses. Detailed modeling of the 
light-curves also provides radii measurements for many of 
these targets. Here, we investigate the properties of the 
EB Kepler-16A & B (|Dovle et al.ll201lL hereafter, Dll), 
which hosts a circumbinary planet Kepler- 16b. Dll 
used a photometric-dynamical model (hereafter, PDM) 
to solve the Kepler light-curve data for the Newtonian 
motion of the three body system, and derived precise 
masses (< 0.5% for the stars, and < 5% for the planet) 
and radii (~ 0.3% for all bodies). Our objective here is to 
use the traditional RV based approach to independently 
measure the masses for Kepler-16A & B. 

Fundamental measurements of a diverse sample of stars 
are critical for testing theoretical stellar mass-radius re- 
lationships, but placing meaningful constraints on mod- 
els requires measur ement precision of better than 2-3% 
(jTorres et all |2010f ). The DEBCa10 catalog fists the 
physical properties of detached EBs with mass and ra- 
dius k nown to better t han ^ 2%, following the crite- 
ria of lAndersenI ()1991[ ) , and currently includes about 
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^ Based on observations obtained with the Hobby- Eberly Tele- 
scope, which is a joint project of the University of Texas at 
Austin, the Pennsylvania State University, Stanford University, 
Ludwig-Maximilians-Universitat Miinchen, and Georg-August- 
Universitat Gottingen. 

^ http:/ /www. astro. keele.ac.uk/jkt/debcat/ 



130 such systems. Most DEBCat stars are solar mass 
or greater, and models and observations agree in this 
regime. Constraints on the low-mass population are 
much sparser: DEBCat includes only 28 stars with 
M < 0.8 Mq, and only 3 with M < 0.3 Mq. Many of 
these have radii 10-20% larger than predicted by inodels 
(e.g.,|L 6pcz-Morales| l2W7l:lL6pez-Morales fc Shawll2007l: 
[Torres et al. 2010). Current theories propose that in- 
creased magnetic fields due to tidal locking produce 
starspots, resulting in decreased effective temperature. 
In addition, the magnetic fields can suppress convec- 
tion in the outer atmosphere of the star. Both processes 
can yield inflated stellar radn rela tive to current models 
(IRibas 2006 HChabrier et al.ll2007D . 

Adequately testing these theories requires expanding 
the sample of low-mass stars with precisely measured 
masses and radii, and spectroscopic observations of the 
Kepler EBs as SB2s can yield such a sample. Because the 
Kepler EBs have a wide range of orbital periods, P, and 
are relatively faint, a comprehensive RV survey is im- 
practical. Various groups are focusing on sub-samples, 
such as transiting hot c ompact objects ( j Coughlin et al.l 
120 111) or low- mass stars (jRowe et al.ll2010( ). We are tar- 
geting ~100 Kepler EBs with light curves that predict 
low-mass secondaries. Our objective is to measure the 
dynamical masses of these stars to 3% or better. Binaries 
with mass ratio, q = M2/M1, much less than one, have 
a secondary-to-primary flux ratio, a, that is very small 
in visible light, but more favorable in the infrared. Such 
systems can often be efficiently solved as SB2s by combin- 
ing visible light spectroscopy that precisely measures the 
primary RVs, with a smaller number of infrared obser- 
vatio ns that detect the secondary (e.g.. iBender fc SimonI 
l2008h . 

We are obse rving our EBs with the H-band APOGEE 
spectrograph (jWilson et al.ll2010l ). which is part of th e 
Sloan Digital Sky Survev-III ([Eiscnstci n et al.l 120111) . 
apogee's multi-object capabilities can survey a large 
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number of EBs as SB2s, but provide a small number 
of epochs for each (typically from three to six). We 
are supplementing these observations with visible light 
spectroscopy from the fi ber- fed High Resolution Spectro- 
graph (hereafter, HRS: [^119981) on the 9.2m Hobby - 
Eberly Telescope (hereafter, HET; iRamsev et all 119981) 
The HET queue-scheduled operation (jShetrone et al.l 
|2007[) . combined with the precise P measured by Ke- 
pler, provides efficient observations of our EB sample at 
targeted orbital phases. In most cases our HRS spectra 
provide single-lined binary orbits (SBls), but some, in- 
cluding the results presented here, have the sensitivity to 
directly solve an EB as an SB2. 

Here, we present new observations of the circumbinary 
planet system Kepler-16 as an SB2. Section 2 describes 
the Kepler-16 system and discusses the importance of cir- 
cumbinary planet systems. Sections 3, 4, and 5, present 
our observations, data analysis, and orbit fitting pro- 
cedures. Section 6 compares our results with Dll and 
draws some relevant conclusions. Future papers will de- 
scribe our ongoing survey in more detail and present dy- 
namical results from additional EBs. 

2. CIRCUMBINARY PLANETS & THE KEPLER-16 SYSTEM 

Planets orbiting binary stars, or circumbinary plan- 
ets, likely provide key insights towards understanding 
planet formation processes, so discovering and charac- 
terizing these systems is a high priority. However, such 
planets are not well suited to discovery with precision 
RV techniques that have yielded most of the confirmed 
planet population because measuring precise velocities 
in a binary spectrum is challenging. Attempts to cir- 
cumvent this problem by combining RVs with interfer- 
ometric observations have resulted in some of the most 
precise s tellar masses known, but have not ye t detected 
planets (|Konackill2005bl: iKonacki et al.ll2010[ ). Measur- 
ing transit timing variations of EBs has become the tech- 
nique of choice in circumbinary planet searches. Eclipse 
timing variation signals, interpreted as due to a plane- 
tary mass object, have been detected around HW Vir- 
ginis, CM Dra, SS Serpens, and HU Aq, although the 
planetary interpretation has also been que stioned (e.g., 
iHorner et alll2011l: Iwittenmver et al.ll2012D . 

Three circumbinary planets have recently been identi- 
fied with Kepler photometry, and characte rized through 
the application of a three-body PDM ([Carter et al.l 
[20T1I): KeDler-16b (Dll), Kepler-34b, and Kepler-35b 
(jWelsh et alll2012[ ). The Kepler-16 system is particu- 
larly interesting for constraining stellar mass-radius re- 
lationships because its stellar components are both low- 
mass. From the PDM, Dll determined that Kepler-16 
has q = 0.2937 ± 0.0006 (J. Carter, private communi- 
cation). When combined with single-lined TRES spec- 
troscopy, they derived masses of Ma = 0.6897 ± 0.0035 
and Mb = 0.20255±0.00065, for Kepler 16-A & B respec- 
tively. Their spectroscopy was unable to detect Kepler- 
16B, which has a = 0.015 in the Kepler bandpass. 

Given the small, but significant, discrepancies between 
stellar radii measured for low-mass stars and current 
stellar models, obtaining robust absolute mass measure- 
ments for individual systems is critical to anchor the ab- 
solute offset for the models. The PDM provides an esti- 
mate of stellar masses based on a three-body Newtonian 
model, and is a relatively new and powerful technique 



that may be applicable to numerous systems with pre- 
cise photometry and multiple transiting objects. An im- 
portant part of scientific investigation is the independent 
confirmation of results and techniques, preferably using 
different instruments and measurements of independent 
variables. Observations of a binary as an SB2 use a Kep- 
lerian model (hereafter, KM) to estimate q and absolute 
masses, and are a well established technique. The SB2 
observations of Kepler-16 we present here directly test 
the PDM results. 

3. OBSERVATIONS WITH THE HOBBY-EBERLY 
TELESCOPE 

We obtained six observations of Kepler-16 with the 
HRS, from September through October 2011, using the 
316g7940 cross-disperser, a 2" diameter fiber that en- 
compassed both stars, and a slit providing R'~30,000. 
Each observation was a 1200 s integration, which yielded 
S/N~300 per resolution element, and was bracketed 
before and after with a ThAr hollow-cathode lamp 
exposure for wavelength calibration. We used the 
same procedure on two known single stars, HD17230 
and GJ905, to serve as spectral templates for Kepler- 
16A & B, respectively. HD17230 has spectral type 
K6V, with Toff = 444 2 ±44K, logq = 4.89 ± 0.06, and 
[M/H] = 0.02 ± 0.03 (IValenti fc Fisheii l2005l) . GJ905 
has spectral type M5 5V, with Teff = 28 00 ± lOOK, and 
[M/H] = 0.00 ± 0.25 (|Leggett et al.l[2000h . 

We extract HRS spectra using a custom opti- 
mal extract ion pipeline, similar to that described by 
ICushing et al.l (j2004j) . written by us in the Interactive 
Data Language (IDL). This pipeline automatically per- 
forms basic image processing, including overscan correc- 
tion, bias subtraction, and fiat fielding. It then auto- 
matically traces the spectral orders in each HRS target 
observation, computes the optimal fiber profile for each 
order, and carries out the extraction. ThAr calibration 
spectra are extracted in a similar manner, using a fiber 
profile determined from exposures of twilight sky or a flat 
lamp observed through the target fiber. A multi-order so- 
lution equating pixel position to wavelen gth is fit to th e 
ThAr spectrum using the linehst of Murph v et all (|2007|) , 
and then applied to the corresponding target spectrum. 
We found no significant difference between the two ThAr 
frames that bracketed each observation, indicating that 
HRS was stable throughout each integration, and so arbi- 
trarily used the later ThAr solution for each observation. 
To facilitate use by other investigators, we have placed 
our extracted Kepler-16 spectra in the NASA Exoplanet 
ArchiveQ. 

We used a synthetic telluric absorption function 

G enerated with t he LBLRTM atmospheric model 
Clo ugh et al.ll2005l) to identify spectral regions contam- 
inated by telluric absorption, and excluded them from 
subsequent analyzes. We removed sharp night-sky OH 
emission features by linearly interpolating over them. In 
the 316g7940 configuration many orders overlap off-blaze 
wavelengths. We combined these orders into five con- 
tiguous segments which were free of telluric absorption: 
6180-6270A, 6340-6450A, 6600-6850A, 7440-7580A, and 
8460-8870A. In each bandpass the K epler-16 binary has 
a different a. Using BT-Settl models (jAllard et al.ll201lD 

^ http:/ /exoplanetarchive. ipac.caltech.edu 
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Table 1 

HRS Radial Velocity Measurements of Kepler-16A & B 



UT Date 


BJD - 2,400,000 


Orbital Phase 


VA(kms-i) 


VB(kms-i) 


2011 Sep 27 


55831.677860 


0.607 


-39.175 ±0.026 


-15.39 ± 0.80 


2011 Oct 6 


55840.643682 


0.825 


-47.615 ±0.025 


13.32 ± 1.03 


2011 Oct 8 


55842.631447 


0.874 


-46.510 ±0.030 


8.88 ± 1.08 


2011 Oct 15 


55849.627832 


0.044 


-30.288 ±0.023 


-44.34 ± 0.63 


2011 Oct 20 


55854.601732 


0.165 


-21.244 ±0.026 


-75.20 ± 0.82 


2011 Oct 24 


55858.588560 


0.262 


-20.917 ±0.025 


-76.80 ± 0.78 



th at correspond to th e Toff, \og{g), and [Fe/H] derived 
bv lWinn et al.l ()2011|) and Dll (and shown in Fig 1) we 
estimate the a in each bandpass to be 0.00522, 0.00743, 
0.00780, 0.0228, and 0.0286, respectively. 

4. SB2 ANALYSIS 

The cross-correlat ion algorithm TODCOR 
(jZucker fc Mazehl I1994D simultaneously cross-correlates 
two template spectra against a target spectrum con- 
taining the blended light from a binary to disentangle 
the component RVs. We used TODCOR, along with 
the HD17230 and GJ905 templates, to measure the 
RVs of Keple r - 16 A &: B in our HRS spectra. Following 
iZucker et al.l ()2003() . we solved each bandpass inde- 
pendently with a locked to the corresponding value 
estimated in [J3l and combined the resulting correlation 
surfa ces with a maximum-likelihood analysis (jZucked 
12003( 1. TODCOR is more sensitive to the RV of each 
component than to the a, and in high-contrast binaries 
where the flux from the primary dominates, a discrep- 
ancy of a factor of 2-3 between the chosen a and the true 
value can still yield equivalent RVs. Similarly, small mis- 
matches in the target and template metallicity primarily 
affect the ability of TODCOR to deriv e an optimal 
a, and do not affect the measured RVs () Bender et al.l 
[2005). By locking a across all HRS epochs, we minimize 
the effect on the RVs of discrepancies with the BT Settl 
a or template metallicity mismatch. 
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Figure 1. BT-Settl models corresponding to Kepler-16A & B, 
plotted along with the Kepler response function. The a of these 
models in the Kepler bandpass, ~ 1%, is consistent with the value 
measured by Dll; in the HRS bandpasses a ranges from ~ 0.5% 
to 3%. 



In analyzing each observation we restricted the sec- 
ondary velocity search to the range bounded by g = 1 and 
the minimum-mass from the mass- function, which are 
constrained by the SBl solution and the orbital phase, 
and chose the largest amplitude peak in this range. The 
component RVs, Va and Vb, are measured by fitting 
the top six to eight points of this peak with a quadratic. 
Uncertainties are derived from the maximum-likelihood 
formalism of Zuckcr (2003), which accounts for spectral 
bandwidth, correlation peak sharpness, and the spectral 
line S/N, but may exclude systematic uncertainties due 
to the instrument or templates. Our templates share 
many of the same spectral features, so the M6V template 
correlates reasonably well with the Kepler- 16 A spectrum 
and results in a spurious peak in secondary velocity equal 
to Va. In all observations except for 2011 Oct 15, the 
true secondary peak was well separated from the spuri- 
ous peak and could be measured directly. For 2011 Oct 
15, the secondary peak was blended in the blue wing of 
the spurious peak. We fit the unblended portion of the 
spurious peak with a Gaussian and subtracted it, which 
provided a clean secondary peak to measure Vb- 

Table [T] lists the midpoint Barycentric Julian Date 
(BJD), the corresponding orbital phase, and our mea- 
sured RVs for each HRS observation. Figure [5] shows cuts 
through the correlation surfaces, along with the mea- 
sured RVs. We are monitoring the HRS long-term sta- 
bility by observing stars known to be intrinsically stable 
to a few ms^^. These observations, taken over many 
months, show RV variations with RMS'^20ms~^, which 
is less than the Va uncertainties reported in Table [TJ 
The secondary peaks for 2011 Oct 06 and 2011 Oct 08 
are weaker than for the other dates (Figure [2]); this is 
likely due to lower intrinsic S/N, and is reflected in the 
uncertainties reported in Table [TJ 

5. DERIVATION OF ORBITAL PARAMETERS AND 
DYNAMICAL STELLAR MASSES 

Measuring Kepler-16 as an SB2 allows us to use a 
KM to derive a dynamical mass-ratio, q, and dynami- 
cal masses, Ma and Mb, for the stellar binary. These 
values constrain the stellar mass-radius relationship and 
test the 23 parameter PDM used by Dll. For reference. 
Table [21 lists the relevant orbital parameters from Dll. 
For our purposes, q is the most fundamental parameter 
describing the EB, and can be deri ved without a ny or- 
bit modehng using the procedure of lWilsoiil ()1941[ ). The 
"Wilson Test" fits a straight line through a plot of Va 
versus Vb: the negative slope of this line directly yields 
q = 0.2986 ±0.0031. 
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Table 2 

Kepler-16 Orbital Parameters and Stellar Masses 



-100 -75 -50 -25 
Primary Velocity (km s"^ 



25 



50 



q 

Q 

O 
O 



0.9655 
0.9650 
0.9645 
0.9668 
0.9664 
0.9660 
0.9508 
0.9504 
0.9500 
0.9708 

0.9702 

0.9696 
0.9710 

0.9705 
0.9700 
0.9728 
0.9724 
0.9720 
-1 
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Figure 2. Cuts through correlation surfaces vs. velocity, showing 
(a) an example primary peak, and (b) secondary peaks for each 
observation. Primary cuts not shown have the same general shape 
as (a). Vertical dashed and solid lines indicate the measured Va 
and Vb, respectively. Dotted lines show the systemic velocity, 7^. 
The spurious secondary peak at Va=Vb (0 ha-s been removed 
for clarity. 

The Kepler photometry timing constrains the orbital 
period, P, and the time of transit, Tt, with precision 
that cannot be improved by new spectroscopy. Table [2] 
lists P and Tt from Dll that we use throughout our 
analysitQ. We expand on the "Wilson Test" results by 
solving Kepler's equations for the HRS RVs, while fixing 
P and Tt to their a priori values, yielding eccentricity, 
e, longitude of periastron, w, semi-major axes, Ka and 
Kb, and systemic velocity 7. We propagated the Dll 
uncertainty on P throughout the analysis, but simplified 
the procedure by using the Dll ephemeris Tt closest to 
our HRS observations and neglecting its uncertainty; the 
impact on the final q, Ma, and Mb should be negligible. 
Dll derived precise values for e, and w, but because these 
are osculating parameters we chose to re-solve for them. 

The TRES RVs reported by Dll constrain the spec- 
troscopic orbit of Kepler-16A. We used the Levenb erg- 
Marquardt fitting code MPFIT (jMarkwardtl I2008D to 
solve Kepler's equations for these RVs, yielding a purely 
spectroscopic set of parameters from the Dll results: 
e = 0.15983 ± 0.00085, w = 263.268° ± 0.040°, Ka = 
13.666 ± 0.012 kms-i, ja = -32.7765 ± 0.0061 kms^^, 
with reduced chi-squared of — 1.11 and 15 degrees-of- 
freedom (DOF). An independent fit to our HRS primary 

The Dll Supporting Online Material incorrectly states that Tt 
is relative to BJD 2455000; actually, the values are given relative 
to BJD 2454900. 



Parameter 



Value 



Photometric-Dynamical Model (Dll) 



P (days) 
Tt 

e 

uj (deg) 
i (deg) 
a (AU) 
9 

Ma (Mq) 
Mb (Mq) 
Ra (Ro) 
Rb (R©) 



41.077580 ±0.000008 
2454965.657623 ± 0.000058 

U.lt)944_o (,Qog2 

263.464too26 
90.340lt°:°°i^ 

n 00/1 Q1 -1-0.00035 

U.ZZ^Oi_Q QgQ34 

0.2937 ±0.0006 
6897+° °°^^ 

U.UO»( _Q QQ34 

n 9nOf;c;+0 00066 

0.6489t°;™}| 

n 90ROQ-I-0.00059 

U.ZZUZO_Q QQQgg 



SB2 Spectroscopy (this work) 



e 

UJ (deg) 
Ka (kms-i) 
Kb (kms-i) 
7^ (kms~^) 
7b (kms^i) 

q 

Ma (Mo) 
Mb (Mq) 



0.15894 ±0.00079 
263.287 ±0.041 
13.642 ± 0.010 
45.56 ± 0.47 
-33.7551± 0.0064 

-33.32 ± 0.35 
0.2994 ± 0.0031 

0.654 ±0.017 
0.1959 ±0.0031 



RVs gives e = 0.1602±0.0020, w = 263.46°±0.12°, Ka = 
13.598±0.022kms-\ and 7^ = -33.770±0.018kms-\ 
with xl = 1-58 for 2 DOF. The HRS and TRES reference 
frames are offset by ~ Ikms^^, which can be accounted 
for with an offset parameter; otherwise, the two solutions 
are in agreement, considering the paucity of HRS RVs. 

We then simultaneously fit the HRS and TRES pri- 
mary RVs and the HRS secondary RVs, accounting for 
the reference frame offset. The number of HRS secondary 
RVs is about four times fewer than the combined num- 
ber of primary RVs, and they have much lower preci- 
sion, so they do not further constrain e and w. Ta- 
ble [2] lists the resulting orbital parameters and mass 
ratio, q = 0.2994 ± 0.0031, which is completely consis- 
tent with q derived from the "Wilson Test" , and Fig [3] 
shows the corresponding SB 2 RV curves. This solution 
has xl = 1.013, with 25 DOF. Our velocity precision 
is sufiicient to require different systemic velocities, ^a 
and 7b, for each star, due to the combined convective 
blue-shift and gravitational redshift of the two stars and 
templates (|Pourbaix et alJ 120021 : [Pravinsl 119991 ) . We es- 
timate the combined effect at ^ 400 ms~^, but did not 
attempt to disentangle the individual contributions. The 
uncertainty measured for 75 implies that 7^ and 75 dif- 
fer by only ~1.25(t. 

To examine the influence of systematic differences be- 
tween the TRES and HRS RVs, we re-computed the or- 
bital parameters using only the HRS data. As before, e 
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and Lo are dominated by the primary RVs. This yields 
e, uj, Ka, and 7^ that are numerically equivalent to the 
HRS SBl derived above, and Kb and 75 that are equiv- 
alent to the values in TableEl with xl = 1-12 for 6 DOF. 
The resuhing q = 0.2986 ± 0.0032 is identical to that 
derived from the "Wilson Test" , and demonstrates that 
nearly all of the reported uncertainty is contributed by 
the secondary RVs: including or excluding the TRES 
RVs does not meaningfully change our derived masses or 
their precision. 

Our measured Ka and Kb yield q ~ Ka/Kb ~ 
0.2994 ± 0.0031 and, using t he ph ysical constants sug- 
gested bv lHarmanec fc Prsal (|2011| ). Ma sin^ i = 0.654 ± 
0.017 Mq and Mssin^i = 0.1959 ± 0.0031 Mq. When 
combined with the inclination measured by Dll, we can 
solve directly for the dynamical masses of Kepler-16: 
Ma = 0.654 ±0.017 M© and Mb = 0.1959 ± 0.0031 Mq. 

6. DISCUSSION 

We have measured q to ~ 1%, consistent with Dll 
at ~ 2tT, and Mb to ~ 1.5%. We performed several 
tests to explore the q discrepancy and the robustness of 
our derived orbital parameters. Using the Dll orbit, 
we calculate predicted Ka — 13.661 ± 0.030 km s^^ and 
Kb = 46.512±0.076kms"^ The difference between this 
Ka and our measured value reported in Table [5] cannot 
account for the q discrepancy. The difference between 
the predicted Kb and our measured value is 2(T. An 
error in our measured Kb of 1 kms~^, consistent with 
our reported uncertainty at 2fT, would account for the 
discrepancy. However, the Ob — Cb residuals (Fig [3]) 
suggest that our Vb uncertainties are overestimated. 

Due to the small number of Vb measurements, a sys- 
tematic error in one might skew our measured Kb- To 
test this, we repeated the KM fit six times, excluding 
each Vb in turn. This produced a range in Kb from 
45.36 ± 0.59 kms~^ to 45.71 ± 0.61kms~^, correspond- 
ing to q from 0.2984 ± 0.0040 to 0.3007 ± 0.0039, and 
7s from -33.21 ± 0.40 to -33.46 ± 0.39. These distri- 
butions are symmetric around the values in Table [21 and 
confirm that no single HRS SB2 observation is skewing 
the measured q. 

We considered that our analysis of the SB2 RVs might 
be affected by the planet's orbital motion. However, this 
is unlikely because Kepler- 16b has a nearly circular orbit 
with P ~229 days (Dll), while our HRS observations 
were obtained over just 28 days. None of our HRS obser- 
vations coincided with a transit event, whic h would sub- 
ject t he RVs to a Rossiter-McLaughlin shift ()Winn et al.l 
120111) . Sunlight scattered off the moon might confuse 
our SB2 analysis, but the moon was down for three ob- 
servations and at least 70° from Kepler-16 for the other 
three. No scattered light was detected in HRS sky fibers 
exposed coincident with our observations, so moonlight 
contamination is unlikely. 

We have used dynamical SB2 measurements to derive 
g, Mai and Mb for the Kepler-16 EB, and have con- 
firmed the q reported by Dll from their PDM to ~2%. 
In addition, the e and w measured by the two techniques 
are equivalent. The claimed precision (g ~ 1% by us 
and ~ 0.2% by Dll), results in a ~2(t disagreement that 
we cannot account for. It is possible that our measure- 
ment errors are underestimated due to an unidentified 
systematic effect, or that the Dll dynamical analysis is 
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Figure 3. RV vs. orbital phase for Kepler-16A & B. Upper Panel: 
Solid circles and squares show our HRS Va and Vb, respectively. 
Open circles show the Dll TRES Va, offset to the HRS RV frame. 
The solid and dashed lines show the best orbital solution (Table 2) . 
Lower Panels: Observed - calculated residuals for the HRS and 
Dll RVs. 

erroneous, potentially due to the presence of additional 
bodies in the system. However, the latter is difficult to 
understand in light of the continued accuracy of transit 
and eclipse ephemerides in subsequent Kepler quarters 
of data (J. Carter, private communication). In addition 
to verifying the PDM result, we have demonstrated that 
masses with 1-2% precision can be measured for Kepler 
EBs, even when no planet transits are present. 
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